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Abstract 

We study the effect of surface fields on the interfacial properties of a 
binary polymer melt confined between two parallel walls. Each wall at- 
tracts a different component of the blend by a non-retarded van der Waals 
potential. An interface which runs parallel to the surfaces is stabilized in 
the center of the film. Using extensive Monte Carlo simulations we study 
the interfacial properties as a function of the film thickness, the strength 
of the surface forces and the lateral size over which the profiles across 
the film are averaged. We find evidence for capillary wave broadening of 
the apparent interfacial profiles. However, the apparent interfacial width 
cannot be described quantitatively by a simple logarithmic dependence 
on the film thickness. The Monte Carlo simulations reveal that the sur- 
face fields give rise to an additional reduction of the intrinsic interfacial 
width and an increase of the effective interfacial tension upon decreasing 
the film thickness. These modifications of the intrinsic interfacial prop- 
erties are confirmed by self-consistent field calculations. Taking account 
of the thickness dependence of the intrinsic interfacial properties and the 
capillary wave broadening, we can describe our simulation results quanti- 
tatively. 

1 Introduction. 

Interface and surface properties in polymeric composites have attracted abiding 
interest. The surface structure of a material correlates with its application prop- 
erties (e.g. adhesive and wetting properties, lubrification) . Hence, the control 
of the surface and interfacial structure is a key for tailoring novel materials. jj| 
Surface and thin film properties also pose challenging theoretical questions. 
The presence of the surface alters the phase behavior^], J| and restricts the 
fluctuations of an interface in the vicinity of the surface, pj However, a direct 
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comparison of surface and interfacial profiles between theory and experiment or 
computer simulation is difficult :[|| While the theoretical approaches calculate an 
intrinsic profile, experimental profiles are broadened by capillary waves. |6|, [?], || 
These fluctuations of the local interfacial position give rise to a dependence of the 
width of the apparent interfacial profiles on the lateral resolution and the thick- 
ness of the film. Though capillary waves are present at all fluid interfaces, |], |lo| 
the study of polymeric systems is particularly rewarding: The extended struc- 
ture of the macromolecules facilitates the application of several experimental 
methods; e.g., the surface excess of a component is typically an order of mag- 
nitude larger than in its low molecular weight counterpart. By virtue of their 
extension a polymer interacts with many neighbors and composition fluctua- 
tions in the bulk are small. Hence, the thermodynamical properties of high 
molecular weight blends are well described by mean field theory and, apart from 
capillary wave fluctuations, self-consistent field theories [|l2[ [ll| |l4|, [ll| are be- 
lieved to provide a detailed description of interfaces in the limit of very long 
chain lengths. Hence, the effects of interfacial fluctuations in polymer blends 
can be accurately isolated and a comparison between experiments, computer 
simulation, and analytic theory is possible. 

Many theoretical approaches j| |) assume short range interactions between 
the monomeric units and the surfaces. Experiments indicate that this is a rea- 
sonable approximation if the range of the interactions is small compared to the 
polymer's extension. |T(| In the case of short range surface fields, the width of 
the apparent profile in a thin film increases like v^D, where D denotes the thick- 
ness of the film. This has been studied analytically, experimentally |l8[ 
and in Monte Carlo simulations. |5|, |l9L In simple fluids, the influence of long 
range forces can shift the order of the wetting transition from second to first 



order 1 21, g2| and modifies the composition profiles near the surface. p3| There 
is also experimental evidence that long range force influence the dynamics of 
wetting layers in thin films, pij Spinodal dewetting in the presence of long range 
forces was recently investigated numerically by Puri and Binder. J25| Effects of 
long range forces have been studied by Monte Carlo simulations by Pereira and 
Wang, |2(| and in the framework of self-consistent field calculations by Genzer 
and co-workers. |27j Both techniques reveal only qualitative changes due to the 
presence of long range surface fields for long polymer chains. The presence of 
long range forces, however, is expected to have pronounced effects on the depen- 
dence of the interfacial width in a thin film on its thickness. In contrast to the 
short range case, a much weaker, logarithmic dependence on the film thickness 
is predicted. This has been recently observed in experiments of Sferrazza et 
al§ 

The present study aims at investigating the interplay between long range 
forces and the apparent interfacial profiles in a thin film geometry. We study 
the dependence of the interfacial width in a binary polymer blend on the film 
thickness and the strength of the surface interactions via Monte Carlo (MC) 
simulations in the framework of a coarse grained lattice model. The results are 
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compared to self-consistent field (SCF) calculations without adjustable parame- 
ter. This allows to separate the effect of capillary fluctuations and the influence 
of the surface interaction on the intrinsic interfacial profile. We restrict ourselves 
to the case of antisymmetric surface fields. I.e., the film is confined between 
two walls; the right wall attracts species A with the same strength and range 
as the left wall attracts species B. Otherwise the two walls are identical and 
the polymer mixture is chosen perfectly symmetric (same chain length, stiffness, 
etc.). 

Analyzing the dependence of the interfacial width on the lateral coarse grain- 
ing size|| [2£| and measuring the spectrum of interfacial fluctuations |2^, ||o], [H| 
we find convincing evidence for a broadening of interfacial profiles due to cap- 
illary waves. The surface fields impart a lateral correlation length £11 on the 
interfacial fluctuations which is determined via the direct measurement of the 
correlation function of composition fluctuations and the spectrum of interfa- 
cial fluctuations. £|| grows proportional to the square of the film width D and, 
consequentially, the apparent interfacial width is expected to exhibit a weak 
logarithmic dependence w ~ In D on the film thickness for large lateral system 
extensions L ^> £m . For the system sizes studied, however, our simulations reveal 
an additional coupling between the long range surface fields and the intrinsic 
properties of the interface. The effect of the surface fields on the intrinsic inter- 
facial properties is confirmed by self-consistent field calculations. A quantitative 
analysis of our Monte Carlo results has to incorporate both effects. 

Our paper is arranged as follows: The next section provides the pertinent 
theoretical background. In section III we describe the model that is used in the 
calculations and briefly comment on the Monte Carlo (MC) and self-consistent 
field (SCF) technique. Then we detail our results and close with a summary 
and a discussion of our findings. 

2 Background. 

In this section, we briefly review some of the basic theoretical predictions and 
provide the analytical description pertinent to the analysis of our Monte Carlo 
results. We consider a binary blend in a thin film of thickness D. The confining 
impenetrable surfaces are parallel to the xy plane and located at z = and z = 
D. Each surface preferentially attracts a different component via a long range 
potential. The temperature is chosen such that the blend is phase separated 
in the bulk (T < T c b) and the components wet their corresponding surfaces 
(T > T wet ) in a semi-infinite slab geometry. In a thin film with antisymmetric 
boundaries, an interface is stabilized in the middle of the film which runs parallel 
to the surfaces. Since the interface is located on average at the center of the 
film the system is in the one phase region (soft-mode phase) . [p"7| 

Let u(x, y) denote the deviations of the local interfacial position from its 
equilibrium position z = D/2. The free energy of a configuration can be de- 
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scribed by an effective interfacial Hamiltonian in terms of the local interfacial 
position:^ 

g4 = /d*d» {?f(Vuf + V cS (u)}*Jd X dy ^f{Vuf + 

(1) 

In the following we measure energies in units of fc^T. The first term accounts 
for the free energy costs of an increase in the interfacial area caused by a fluc- 
tuating local interfacial position. For an unconfined interface a e g is expected 
to equal the interfacial tension a between coexisting bulk phases. This has 
been confirmed by Monte Carlo simulations for the model employed in the 
calculations. (3(| The second term describes the effect of the surface fields on 
the local interfacial position. For small deviations u the potential can be ap- 
proximated by a parabola. Fourier decomposing the local interfacial position 
and using the equipartition theorem, we obtain for the mean square value of the 
Fourier component u(g):H 
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where the parallel correlation length for interfacial fluctuations is given by: 



(2) 



e " V d*v cS {u)/du* (3) 

For surface fields which decay like z~~ n the second derivative of the effective 
potential V c s in the middle of the film is of the order D~ n ~ x (see below). Con- 
sequently, the parallel correlation length increases proportional to £)+("+ 1 )/ 2 . 
This length scale acts as a cut-off for the spectrum of interfacial fluctuations 
and limits the broadening of the apparent profile. The fluctuations of the local 
interfacial position u(x,y) are Gaussian distributed and its variance is obtained 
by integrating over all lateral Fourier components. The power law spectrum 
of the free interface (V e g = 0) leads to logarithmic divergences for q — ► and 
q — > co. We can remove these divergences by introducing heuristic cut-offs g max 
and q n 



(^y)) = ^J* -^^t (4) 




In a system of finite lateral extension L the lower cut-off is simply determined 
by the system size g m i n = 2-k/L. The value of the upper cut-off can be defined 
as the lateral size Bq on which the width of the apparent interfacial profile 
coincides with the width of the intrinsic profile calculated in the self-consistent 
field calculations. |Eq| 
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We approximate the laterally averaged apparent profile by the convolution 
of the intrinsic profile po(z) with the Gaussian distribution of the local inter- 
facial position u(x,y). This approximation is appropriate for an unconfincd 
interface, but it neglects a possible dependence of the intrinsic profile on the 
local interfacial position, i.e., the distance between the interface and the wall. 
Defining the width of the interfacial profile by the inverse slope at the center, 
this convolution approximation yields:^, [l(], 33, 33] 



w = w n + 



dz 3 



1 

8^ 



In 



where we have approximated the shape of the intrinsic profile by an error func- 
tion in the last step in order to calculate the numerical prefactor. 

In the limit -Bo £ii ~ D 2 <C L the apparent interfacial width increases 
like w 2 ~ (l/2cr)lnD with the film thickness for non-retarded van der Waals 
interactions (n = 3). Sferrazza et al. found experimental evidence for such a 
weak logarithmic dependence of the apparent width on the film thickness in 
PS/PMMA polymer blends. |8| However the interfacial tension a extracted from 
the thickness dependence was about a factor 1.6 larger than estimated by other 
means. The aim of our Monte Carlo simulations is to investigate the influence 
of long range surface fields on the interfacial width taking advantage of the 
simultaneous access to the interfacial structure on various lateral length scales 
in the Monte Carlo simulations. 



3 Model and computational techniques. 

Coarse grained models of polymeric systems have proven useful for investigat- 
ing universal properties of polymer blends and solutions. We employ the bond 
fluctuation model |3|| on a three dimensional simple cubic lattice. This model re- 
tains the relevant features of polymeric materials - connectivity of the monomers 
along a chain, excluded volume interaction of the segments, and a short range 
thermal interaction potential between monomer units - and is highly computa- 
tionally efficient. In the framework of the bond fluctuation model each segment 
occupies all eight corners of a unit cell, and no site can be doubly occupied. 
Adjacent monomers along a chain molecule are connected via one of 108 bond 
vectors of length 2, y/E, 3, or y/lO. Here and in the following all lengths are 
measured in units of the lattice spacing. We work at a filling fraction 8$ = 0.5 
which corresponds to a melt or a concentrated solution. ]3q] The binary blend 
consists of two species denoted A and B. Both species are structurally sym- 
metric and comprise N — 32 monomeric units. Since each monomer in our 
coarse grained model corresponds to 3-5 chemical repeat units the chain length 
corresponds to a degree of polymerization of the order 100. 

The monomers interact via a square well potential which is extended over all 
54 neighbors sites up to a distance y/6. This includes all lattice sites inside the 
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first peak of the density-density pair correlations function. A pair of monomers 
of the same type in the range of the square well potential lowers the energy 
by e, whereas unlike species increase the energy by the same amount. We 
choose the well depth e = 0.03, a value which corresponds to a rather strong 
segregation \N ~ 5.1 where \ denotes the Flory Huggins parameter. p7[ Many 
properties of this coarse grained model have been determined in prior studies: 
The radius of gyration is R g — 6.96. The interfacial tension a = 0.015 has been 
measured via reweighting techniques Q and it gratifyingly agrees with the value 
measured from the spectrum of interfacial fluctuations of a free interface. |3(J The 
correlation length of composition fluctuations |5) in the bulk is 3.8 and the 
intrinsic width of the free interface is 4.6 in self-consistent field calculations. []30|, 

n 

The confining surfaces are modelled as ideally flat, parallel, and structureless. 
The left wall attracts the A component whereas the right one the B species. 
The long range potential between the monomers and the wall is taken to be 
i£iuVb(z), where the positive sign holds for A monomers and the negative for 
B ones. The shape of the potential is given by 

1 for z < A 

^ (; ) { i - (D-2-zy lox \<z<D-2-\ (6) 
-1 for D-2-\<z 

as to mimic non-retarded van der Waals interactions between the monomers 
and the surface molecules. The parameter combination e„,A 3 corresponds to 
the Hamaker constant. The monomer coordinates refer to the left corner of the 
monomers, hence the z coordinates range of a monomer range in the interval 
[0,D-2]. The left wall attracts B monomers at position z as much as the right 
wall attracts A ones at position D — 2 — z, i.e., the surface fields are completely 
antisymmetric. This potential shape is illustrated in Fig.|l|. The strength of the 
potential is held constant at e w = 0.1, whereas we vary the range 1 < A < 15. 
Prior studies]5|, [u| show that at this combination of parameters the A phase 
wets the surface in a semi-infinite geometry for purely short range surface fields 
which are extended over the two nearest layers to the wall. Since the long range 
surface fields are always stronger, the system is above the wetting temperatures 
for all values of A. In a thin film geometry with antisymmetric surface fields an 
interface is stabilized which runs parallel to the walls and is located on average 
in the center of the film. jl7j 

Simulations are performed in the canonical ensemble and in the semi-grandcanonical 
ensemble. To relax the polymer conformations on the lattice we use a combi- 
nation of random local monomer displacements and slithering snake moves. |j37| 
While the composition of the system is fixed at pa = 0.5 in the canonical runs, 
the composition fluctuates in the semi-grandcanonical simulations and the ex- 
change chemical potential Afj, between the polymer species is controlled. By 
virtue of the symmetry with respect to exchanging the monomer labels A ^ B 
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and the direction of the z axis, the average interfacial position fluctuates around 
the center of the film at A/j, = 0. Though all other than local hopping moves do 
not correspond to a physically realistic dynamics they yield a reasonably fast 
equilibration. The different Monte Carlo moves are applied in the ratio: local 
hopping : slithering snake : semi-grandcanonical moves = 4:12:0.4 as to equili- 
brate the polymer conformations and the composition fluctuations on roughly 
the same time scale. 

To investigate the interfacial structure on different lateral length scales (block 
analysis),)^] we divide the system laterally (i.e. parallel to the interface) into 
subsystems of size B x B. In each of these B x B x D columns we determine 
the local interfacial position u(x, y) using the concept of the Gibbs dividing 
surface. || Then we average the profiles with respect to this position and deter- 
mine its width. Employing this local position of the interface we calculate the 
spectrum of interfacial fluctuations. [p9| [30f 

In order to distinguish between the influence of the long range surface fields 
on the capillary wave broadening and the intrinsic interfacial properties we in- 
vestigate the same model in the framework of the self-consistent (SCF) field 
theory. A detailed account of our self-consistent field (SCF) method has been 
given elsewhere 1 30, 3l[] and we only summarize the salient features pertinent to 



the present study. In the mean field approximation the mteraiolecular interac- 
tions are replaced by potentials of mean field. The potential wa acting on the 
zth monomer of A polymer a at position r a ^ takes the form 

WA(r a ,i) = -ze(l + t^o^) [pA{r a ,i) - PB(r a ,i)] + C [pA(r a ,i) + PB(r a ,i) - 1] 

(7) 

and a similar expression holds for wb- The first term describes the intermolecu- 
lar interaction of an A monomer via the square well potential of depth e with its 
surrounding, z = 2.65 denotes the number of monomers of other chains in the 
interaction range. |37j It has been extracted from MC simulations of the inter- 
molccular paircorrelation function of a dense melt. The derivative in z direction 
takes account of the finite range l Q = 16/9|lJ |§ of the square well potential. 
The expression mimics the "missing neighbor" effect at the surface. |l2| The 
second term represents the free energy cost of deviations of the local monomer 
density pa + Pb from the reference value 1. This free volume contribution is ap- 
proximated by an expression introduced by Helfand, [fl2| and the compressibility 
£ = 4.ljllj has been determined from Monte Carlo simulations of the athermal 
(e = 0) system. 

The local fields depend on the local densities which, in turn, are calculated 
by the average of single chain conformations in the potential of mean field and 
the surface field: 

Ea=i Eilfi V6(r - r a ,i) exp (- J2?=i w A (r 
PA(r) = pa 7 ^— * s '- (8) 

Ea=l^»M eX P (-T,i=l w A(r a ,i) . 
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V w is the Boltzmann weight of an isolated chain with respect to the surface 
fields, i.e., V w vanishes if the z coordinate of a segment is located outside 
the interval [0, D ~ 2] and otherwise takes the form exp(^j = i e w Vo{Zi)). The 
summation runs over a representative sample of up to 15 million single chain 
conformations which have been extracted from bulk simulations of the bond 
fluctuation model at the same density. We evaluate this average on a CRAY 
T3E parallel computer, assigning a subset of conformations to each proces- 
sor. The calculations incorporate the detailed chain structure on all length 
scales without any adjustable parameter. The above equations are expanded 
into a Fourier series |Q and solved self-consistently. We compare our simula- 
tional results quantitatively to the results of the self-consistent field calculations. 
The self-consistent field scheme describes the detailed structure of the uncon- 
fined interface quantitatively, |]30| except for the broadening due to capillary 
waves. It has also been successfully applied to studying wetting and capillary 
condensation. |3l| 

4 Results. 

We illustrate the influence of the surface fields and the film thickness on the 
interfacial profiles in Fig.^(a). It displays the apparent profile of the order 
parameter m(z) = [pa(z) ~ Pb(z)]/[pa(z) + pb(z)], i.e., the profile is averaged 
over the complete lateral extension L = 256 of the simulation cell. Increasing 
the strength of the surface fields A at constant film thickness results in a decrease 
of the width of the apparent interfacial profile for lateral extension L — 256 by 
a factor of 2. The inset presents the results of the SCF calculations for the 
same parameters. The interfacial profiles in the SCF calculations are narrower, 
indicating the broadening of the apparent profiles by capillary waves. The SCF 
data also exhibit a dependence on the potential parameter A which is, however, 
weaker than in the MC simulations. Fitting the order parameter profile to 
the expression m(z) = tanh[(z — D/2)/w], we determine its width w. The 
dependence of the apparent interfacial width on the film thickness D and the 
potential shape A is presented in Fig.||(b). Increasing A we reduce the interfacial 
width. A similar, though weaker dependence is found in the SCF calculations. 
Moreover the apparent width in the MC simulations increases stronger than in 
the SCF calculations upon increasing the film thickness at small and moderate 
values of A. This additional dependence in the simulations can be attributed to 
capillary wave broadening: the larger the film thickness the larger is the lateral 
correlation length and the stronger is broadening due to fluctuations of the local 
interfacial position. 

To proceed in quantifying the effect of capillary waves, we investigate the 
dependence of the interfacial width on the lateral coarse graining size,|^| which 
is set by the length B of the columns used in the averaging procedure. The 
MC results for A = 5 are shown in Fig.||(a) for several film thicknesses D. For 
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the larger film thickness, the interfacial width increases logarithmically with the 
block size. Consequentially, the upper cut-off for the capillary wave broadening 
is set by the lateral coarse graining size <7 m i n = 2ir/B rather than by the thickness 
dependent parallel correlation length £|| . Upon decreasing the film thickness, 
there are deviations from the behavior of a free interface; the interfacial width at 
large block size B tends to a limiting value which depends on the film thickness. 
This shows the limiting of the broadening due to the parallel correlation length, 
which is smaller than the lateral block size B. For film thickness which are 
comparable to the polymer's end-to-end distance there is almost no influence of 
the block size on the interfacial width. 

The vertical line in Fig.||(a) marks the lateral length scale Bq = 8 on which 
the apparent interfacial width of a free interface (5) coincides approximately with 
the results of the SCF calculation wscf = 4.6. This yields an estimate for the 
short wavelength cut-off, beyond which fluctuations are describable via the in- 
trinsic profile. However, even on this length scale the interfacial width decreases 
upon narrowing the film thickness. An effect of similar magnitude is also ob- 
served in the SCF calculations (cf. Fig.^(b)). Hence the surface fields modify the 
intrinsic properties of the interface. A possible dependence of the short wave- 
length cut-off on the film thickness appears to be only a minor effect because of 
the agreement between MC and SCF results. 

In Fig.^(b) we complement this real-space analysis of the apparent interfacial 
width with the measurement of the Fourier spectrum of interfacial fluctuations. 12! 

In accord with the theoretical expectation the inverse mean square amplitude 
of the Fourier components of the local interfacial position depends quadratically 
on the lateral wavevector q. The slope of the curve yields the interfacial tension 
er c ff of the capillary wave Hamiltonian,|30], |3lJ while the abscissa of the line gives 
an estimate for the inverse of the lateral correlation length (cf. eqn.(|^)). Upon 
increasing the film thickness we increase the lateral correlation length and the 
effective interfacial tension approaches the interfacial tension of a free interface 
from above. 

The results of the intrinsic properties for the interfacial width wJBq — 8, D) 
and the effective interfacial tension a e g(D) are summarized in Fig.^. Reducing 
the film thickness to about 2R g , we decrease the intrinsic interfacial width by 
more than a factor 2. This result is in qualitative agreement with the SCF 
calculations. The small quantitative discrepancies between the MC result at 
fixed block size and the SCF calculations are probably due to residual effects of 
capillary waves and difficulties of the SCF calculations to capture quantitatively 
the detailed structure (e.g., packing effects) of the polymeric fluid in the vicinity 
of the surface. J3l| With the decrease of the intrinsic interfacial width goes along 
a pronounced increase of the effective interfacial tension by more than a factor 
2.5. Increasing the potential parameter A has similar effects on the intrinsic 
properties as reducing the film thickness D. 

The thickness dependence of the lateral correlation length is presented in 
Fig.Ba). It can be rationalized as follows. The energy of the interface placed 
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into the long range surface field can be described by an effective potential 
Vlk(u). Assuming that the interfacial profile of the interface centered at z = 
-D/2 + u takes the form pa (z) = 1 — ps{z) = (1 + tanh[u/u>o]/2, we can calculate 
the long range contribution to the energy: 



Vlr(u) w / dz $e w V Q {z) [p A {z) - p B (z)} 



d 2 V hR (u) 192<^e w X 3 ( 10^ 2 / Wo \ 2 , 1Kfl1 fw 



(9) 



where $ denotes the monomer number density. The leading term represents the 
energy of a kink-like interface in the surface field. To this leading order, the 
parallel correlation length takes the form 



^^192^-^ (10) 

As mentioned above, the correlation length increases with the square of the film 
thickness D. 

In addition to the corrections due to the finite width wq of the intrinsic 
profile, the surface fields alter the composition profile in the vicinity of the 
wall. This distortion decays exponentially with the distance from the wall and 
gives rise to a short range interaction between the surface and the interface. 
The effect of this short range contribution can be satisfactorily described by a 
parallel correlation length £|| SR of the form(D 

£||sr « « exp I I (11) 

where £b denotes the correlation length of composition fluctuations in the bulk 
and KT 1 — £(,(1 + l/87r£ 2 cr) the range of the effective interaction between the 
surface and the interface. In the absence of long range surface fields the parallel 
correlation length increases exponentially with the film thickness D. In the 
present case we assume that the effective interactions of the surfaces with the 
interface are additive and the inverse squares of both parallel correlation lengths 
add, respectively: 

«r a ( i2 ) 

In Fig.|| (a) we compare this approximation with the MC results and SCF calcu- 
lations. In the simulations the parallel correlations length has been determined 
via the spectrum of interfacial fluctuations and the direct measurement of the 
composition correlation functions |B 

_ (m(r, z)m(0, z)) - (m(0, z)) 2 
9[r ' Z> - (m(0,z) 2 ) - (m(0,z)) 2 ( ' 
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In the middle of the film we can estimate £|| viag(r,z = D/2) ~ l/y^rexp(— r/£ii). 
The results of both independent measurements agree quantitatively for small D, 
whereas the differences at larger D are due to statistical errors. The figure shows 
that the long range forces alone yield a too large estimate for the parallel corre- 
lation length £|| , whereas the combination of long range and short range surface 
fields gives a reasonable quantitative description of the simulation data. This is 
also confirmed by detailed SCF calculations. In the SCF framework we measure 
the free energy of a film as a function of its composition p = 1/2 + u/D around 
/3=l/2 and employ the bulk value of the interfacial tension a to compute the 
lateral correlation length according to eqn.(|^). For comparison we also present 
the lateral correlation length in the case of purely short range surface. For the 
parameter combination investigated, £sr is always much larger than its long 
range counterpart - a fact which shows that the behavior is dominated by the 
long range surface fields. 

In Fig.||(b) we display the dependence of the parallel correlation length on 
the surface field A for thickness D = 32 and 48. For small surface fields Awl, 
the correlation length approaches the value expected for short range forces, while 
for larger strength of the surface fields the long range contribution dominates, 
and we expect £|| ~ A~ 3 / 2 . The limiting power law behavior is, however, only 
reached for fairly large surface fields (cf. inset). The analytic approximation 
above describes the crossover between the short range dominated (A small) to 
the long range dominated (A large) reasonable well and the behavior is also 
reproduced by our SCF calculations. The SCF calculations agree very well for 
large and moderate values of A, whereas there are deviations for small A. In the 
latter case the short range contribution dominates and the lateral correlation 
length stems from the truncation of the interfacial profile at the wall. The 
deviations between the MC results and the SCF calculations might indicate 
that the SCF calculations capture only qualitatively the fine details of the profile 
at the wall. We present in Fig.|| the ratio between the purely long range (cf. 
eqn.([io|)) and short range contributions (cf. eqn.(jll])). This illustration shows 
that the long range forces are always dominant for large film thickness D whereas 
the regime where short range surface fields dominate the behavior is restricted 
to rather small values of A and film thickness D. 

Having determined the influence of the surface fields on the intrinsic inter- 
facial profiles, we turn to the apparent interfacial width w which is accessible 
in experiments. J|, H, |?| FigJ^(a) presents the dependence of the apparent inter- 
facial width on the film thickness for A = 5. Both data for canonical simula- 
tions and for semi-grandcanonical ones are displayed. Similar to prior studies 
of short range surface fields,^ the apparent interfacial width in the canonical 
ensemble saturates with growing film thickness at a value of the free interface 
broadened by capillary waves. The limiting broadening is determined by the 
cut-off g m i n = 2ir / L rather than the lateral correlation length. In the semi- 
grandcanonical ensemble, the average composition fluctuates, hence the inter- 
face in a thick film of small lateral extension is only weakly bound to the center 
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of the film and undergoes a diffusive motion. This leads to a linear increase of 
the apparent interfacial width with the film thickness. |{| 

Similar to the experimental procedure, we have tried to estimate the intrinsic 
properties of the interface by describing the MC data for the apparent width via 
eqn.(||). UsingQ B = 12 and £y according to eqn.(|Io|), we fit the data (solid 
line) and obtain for the intrinsic width wq = 2.1 and the interfacial tension 
a = 0.022. The curve does not describe the MC results well and the extracted 
values for the intrinsic width and the interfacial tension differs from the SCF 
results by a factor 0.5 and 1.5, respectively. In a second attempt we assume 
£|l = const!? 2 and leave the proportionality constant as a fit parameter. In 
the experimental situation this procedure corresponds to an unknown Hamaker 
constant. The dashed line displays the resulting fit which captures the behavior 
at small film thicknesses, but describes the thick film data less satisfactorily. 
The extracted parameters for the intrinsic width and the interfacial tension 
deviate from the SCF values by a factor 0.3 and 0.7, respectively. Moreover, 
the value of the parallel correlation length is overestimated by more than an 
order of magnitude. Therefore we conclude that our MC results cannot be 
described appropriately by a simple effective interfacial Hamiltonian and it is 
not possible to extract reliable values for the intrinsic interfacial profiles from 
such an analysis. This result emphasizes the importance of the modification 
of the intrinsic interfacial properties by the presence of the long range surface 
fields. 

However, using the measured values for the intrinsic interfacial width w(Bq = 
8,D) extracted via the block analysis, the effective interfacial tension <7 c ff(-D) 
as determined via the capillary fluctuation spectrum, and the parallel correla- 
tion length £\\(D) according to eqn.(p"2|) we can describe the MC results quan- 
titatively with the theoretical prediction (||). This is displayed in Figj^(b) 
where we present the increase of the interfacial width due to capillary waves 
wl ap = w 2 -w 2 (D) = l/8<7 cff (£>) H[27r/B } 2 +q 2 )/({2ir/L} 2 +q 2 ). The long 
dashed line marks the simple logarithmic dependence (l/2<r) In I? for compari- 
son. 

5 Discussion and outlook. 

We have presented detailed Monte Carlo simulations targeted to investigating 
the effect of long range surface fields on the interfacial profiles in thin binary 
polymer films. The effect of capillary waves has been demonstrated conclusively 
via the dependence of the interfacial width of the apparent profiles on the lateral 
coarse graining size B and a measurement of the Fourier spectrum of interfa- 
cial fluctuations. Both the strength of the surface fields A and the thickness 
of the film 2.3R g < D < !8.5R g influence the apparent interfacial width. The 
surface fields give rise to a lateral correlation length £|| which acts as a long 
wavelength cut-off for the spectrum of capillary waves. This parallel correlation 
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length incorporates two contributions - a short range one which stems from the 
distortion of the interfacial profile near the wall and a part caused by the long 
range surface field. For large film thicknesses the long range contribution domi- 
nates and we find £|| ~ D 2 . The short range distortion contributes notably only 
for small surface fields and film thickness (D w 5R g ). This is confirmed by two 
independent measurements in the Monte Carlo simulation and a quantitative 
comparison to self-consistent field calculations. 

However, our Monte Carlo simulations are not compatible with a simple log- 
arithmic dependence of the width on the film thickness, which has been assumed 
in experiments by Sferrazza et aZ.Q Instead, the Monte Carlo simulations show 
an additional influence of the surface fields on the intrinsic interfacial proper- 
ties: upon decreasing the film width to 2.3i? s , the intrinsic interfacial profile 
becomes narrower and the interfacial tension increases by more than a factor of 
2. An effect of the long range surface fields on the intrinsic properties persists 
up to film thicknesses which exceed the radius of gyration by roughly a factor 
of 10. Similar effects have been observed for short range interactions. ||, [nj The 
modification of the intrinsic interfacial properties is confirmed by self-consistent 
field calculations, which yield reasonable quantitative agreement without any 
adjustable parameter. The corrections to the effective interfacial Hamiltonian 
are too strong in our MC simulations to extract the intrinsic interfacial prop- 
erties from a simple logarithmic growth of the apparent interfacial width with 
the film thickness. Taking due account of the capillary wave broadening and 
the influence of the surface interactions on the intrinsic profile simultaneously 
we can describe our Monte Carlo results quantitatively. 

This observation might also offer an explanation why the value of the interfa- 
cial tension extracted from the thickness dependence of thin film experiments by 
experiments of Sferrazza et al. B is larger than estimated via mean field theory. 
Deviations from the behavior described by the effective interfacial Hamiltonian 
have also been found in experiments of Kerle et al. [Q, ^] In both experiments 
the simple theory overestimates the apparent interfacial width if the intrinsic 
(i.e., self-consistent field) results for the interfacial widths and tension are em- 
ployed. Moreover, in these two experiments the range of film thickness are of 
the order of several R g , a situation comparable to the systems studied here via 
Monte Carlo simulations. We believe that a quantitative understanding of the 
interfacial properties in thin films calls for a joint analysis of the interplay be- 
tween surface forces and capillary waves as well as of the influence of the surface 
forces on the intrinsic interfacial properties. 
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Figure 1: Illustration of the long range surface field Vq{z) acting on an A monomer at position z. The 
parameter A tunes the shape of the surface interactions (1 < A < 10). The shaded region corresponds to 
purely short range surface interactions. The thick curves marks the potential for A = 5. 




Figure 2: (a) Apparent interfacial profile in the Monte Carlo simulations (L = 256) and intrinsic profiles 
in the self-consistent field calculations (inset) for D = 32, e w = 0.1, and e = 0.03. The different curves 
refer to A = 1, 5, and 10. (b) Dependence of the apparent interfacial width on the shape parameter A for 
D = 32,48 and 64 and lateral system size L = 256 in the Monte Carlo simulations. The inset displays 
the result of the self-consistent field calculations for the same parameters. 
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Figure 3: Influence of the lateral coarse graining size on the interfacial properties: (a) Dependence of 
the interfacial width on the lateral block size B for A = 5 and different film thicknesses D. (b) Fourier 
spectrum of interfacial fluctuations for A = 5 and different film thicknesses D. 



5 ■ 
4 - 
3 - 
2 - 
1 



32 




64 
D 



96 



128 



Figure 4: Dependence of the interfacial width at lateral block size B — 8 and the effective interfacial 
tension as extracted from the interfacial fluctuation spectrum on the film thickness D at A = 5. The 
open symbols present the result for the interfacial width in the self-consistent field calculations. 
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Figure 5: 

(a) Parallel correlation length as a function of the film thickness D for A = 5. Filled squares correspond to 
the measurement via the correlation function g(r, z = D/2), whereas open circles present the analysis of 
the interfacial fluctuation spectrum. Stars denote the results of the self-consistent field calculations. The 
Monte Carlo results for purely short range surface fields^) are depicted as filled squares. The dashed line 
corresponds to purely long range surface fields, the solid one presents the results of combined short range 
and long range effects according to eqns ( |Io| ) and (|l2|). (b) Dependence of the parallel correlation length 
on the potential parameter A for D = 32 and 48. Filled symbols correspond to Monte Carlo data, whereas 
open ones denote the results of the self-consistent field calculations. Dashed and solid curves correspond 
to purely long range forces and the combination of short range and long range forces, respectively. The 
arrows on the left hand side indicate the values for purely short range surface fields. The inset displays 
the Monte Carlo data on a log-log scale. 




Figure 6: 

Ratio of the correlation length for short range surface fields and long range surface fields according to 
eqns.(^fj) and ([l(]) as a function of the film thickness and the strength of the surface fields A. For large 
film thicknesses D the long range forces always dominate, whereas the effect of short range forces becomes 
appreciable for small D w 5R g and small values of A. 




Figure 7: 

Dependence of the apparent interfacial width on the film thickness D at A = 5. (a) Filled symbols denote 
the apparent interfacial width in the canonical simulations, and the arrows on the right hand side mark 
the width of an unconfined interface at the lateral system size L. Open circles display the results of the 
semi-grandcanonical simulations. The lines correspond to fits without considering the the modification 
of the intrinsic interfacial properties, (b) Broadening w 2 ap = w 2 (L = 256) — w 2 (B = 8) of the interfacial 
width upon increasing the film thickness. The symbols correspond to the same Monte Carlo results as in 
(a), whereas the solid curve presents a detailed description according to eqns. d) and @. The dashed 
line corresponds to a simple logarithmic dependence. 



